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O^, : Abstract 

O I It was argued by Brigante et.al [1] that the lower bound on the ratio of the shear viscos- 

ity to the entropy density in strongly coupled plasma is translated into microcausality 
violation in the dual gravitational description. Since transport properties of the system 
/\f ■ characterize its infrared dynamics, while the causality of the theory is determined by its 

5^ \ ultraviolet behavior, the viscosity bound/microcausality link should not be applicable 

to theories that undergo low temperature phase transitions. We present an explicit 
model of AdS/CFT correspondence that confirms this fact. 
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1 Introduction 

Working in the framework of the gauge theory/string theory correspondence of Malda- 
cena [2,3], Policastro, Son and Starinets computed the ratio of the shear viscosity rj to 
entropy density s of the A/" = 4 SU{N) supersymmetric Yang-Mills (SYM) plasma, in 
the planar ('t Hooft) limit and for infinitely large 't Hooft coupling^ A = Qyu^ ~^ ^^ 
[4], finding 

^ = -. (1.1) 

s An 

Shortly afterwards it was argued in [5] that (1.1) is in fact a universal result in all 
gauge theory plasma at infinite coupling that allow for a dual holographic description^. 
The holographic result (1.1) is remarkable in a sense that a simple quasi-particle 
picture of hydrodynamic transport suggests a quantum mechanical bound [12] 

->Oil). (1.2) 

s 



^We set h^ks = 1. 

^Further generalizations/proofs of the shear viscosity universahty theorem appeared in [6-11] 



This fact, along with the observation that all known fluids in nature^ have larger shear 
viscosity to entropy density ratios, led Kovtun, Son and Starinets (KSS) to conjecture 
a bound for any fluid [6]: 

->^- (1-3) 

It is possible to construct a phenomenological counterexample in which, by increasing 
the number of species in the fluid while keeping the dynamics essentially independent of 
the species type, the bound can be violated [14]. Unfortunately, the particular example 
[14] does not have a well-defined relativistic quantum field theory completion [15]. 

The first test confirming the KSS bound, at least for A/" = 4 SYM at large (but 
finite) 't Hooft coupling, was done in [16]^. The finite 't Hooft coupling corrections 
on the gauge theory side translate into higher-derivative gravitational corrections on 
the string theory side of the holographic correspondence [3] . This model, along with 
generalizations [19,20], describes a superconformal gauge theory plasma — a consistent 
relativistic quantum field theory — with the same anomaly coefficients (central charges) 
c = a in the trace of the stress-energy tensor, 

(T".)c-r=j^h~E,. (1.4) 

Here E4 and I4 correspond, respectively, to the four-dimensional Euler density and the 
square of the Weyl curvature: 

E4 = R.upxR^"'' - ^R^uR"" + R\ h = R^upxR^"'' - 2R,,R^''' + ^R^ . (1.5) 

In [22] Kats and Petrov put forth the first consistent example of a relativistic quantum 
field theory which violates the KSS viscosity bound^ — the M = 2 Sp{N) supercon- 
formal gauge theory plasma with 4 fundamental and 1 antisymmetric hypermultiplets. 
The violation of the viscosity bound can be traced back to the inequality between 
the central charges of the theory, c ^ a. More precisely, the bound is violated once 
c — a > 0, which is generic in superconformal gauge theories with c ^ a [24]. Moreover, 
since c — a ^ N, this is a finite A^ correction, and is not due to having finite 't Hooft 
coupling. Once again, the inequality for the central charges on the gauge theory side 
translates into particular higher- derivative corrections to the supergravity approxima- 
tion [25], which, to insure reliable computations, have to be regarded as being 'small'. 



■^A strongly coupled Quark-Gluon Plasma might be a counterexample [13]. 
^For further analysis and generalizations see [17-21]. 



^See also [23] 



As a result, the KSS bound violation in holographic models realized in string theory is 
necessarily perturbative. 

The work [22] convincingly established that the original KSS bound (1.3) can not 
be a quantitative formulation of a loose quantum- mechanical bound (1.2). Thus, the 
question remained as to whether or not a bound of the type (1.2) existed. As we 
mentioned above, because of the universality of the shear viscosity in the supergravity 
approximation, any finite violation of the KSS bound has to be studied in a holographic 
model of the AdS/CFT correspondence, rather than a particular realization of the 
holographic correspondence in string theory. A simple enough model to fulfill this 
purpose is that of Gauss-Bonnet gravity with a negative cosmological constant [23] : 

1 

Up to field redefinitions, for Xgb ^ 1 the gravitational model (1.6) is equivalent to 
the string theory holographic example of Kats and Petrov [22], for sufficiently large 't 
Hooft coupling, where one identifies 

^ = AXoB + O {Xl,) . (1.7) 

The advantage of (1.6) compared to [22] is that the former gravitational model is 
consistent for arbitrary values of Xgb > \ [23]. As such, it defines via the AdS/CFT 
correspondence a dual conformal gauge theory plasma, which we call GB plasma, with 
central charges [26] 



Sgb = H73" / d^xy/^ 



12 A 



;i.6) 



and hence 
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23/2 £3^ 
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1 + Vl - 4Agb)'/' Vl - 4A, 



GB 



AXan?'^ 3A/l-4Ar7R-2 , (1. 



^ -1). (1.9) 



c \^1-AXgb 

Notice the parallel with the construction of [14]: we identified a relativistic quantum 
field theory as a holographic dual to (1.6), with a shear viscosity to entropy density 
ratio given by [23] 

'^ = ^{1-AXgb). (1.10) 

s Ait 

which apparently leads to an arbitrary violation of the KSS bound (or any bound of 
the type (1.2)), given appropriate choices of Xgb (or equivalently of the central charges 



of the theory). To complete the analysis one needs to address the question of the 
consistency of the GB plasma, as a relativistic quantum field theory. This was done 
in [1,27]. It was found that once 

Ac.>4, (l.U) 

the spectrum of excitations in the GB plasma contains modes that propagate faster 
than the speed of light [1]. Likewise, for 

Agb<-t^, (1.12) 

the GB plasma also contains microcausality violating excitations [27]. Given (1.11) 
and (1.12) we are led to conclude that consistency of the GB plasma as a relativistic 
QFT constrains its viscosity ratio to be® 

16 r? 16 , , 

- > 4.f > -, (03) 

Exactly the same constraint arises by requiring "positivity of energy" measured by a 
detector in the GB plasma [31]. 

To summarize, the example of the GB plasma appears to suggest a link between 
the violation of the shear viscosity bound of the type (1.2) and the violation of micro- 
causality /positivity of energy in the theory^. In this paper we argue that such a link 
can not be of fundamental nature. 

Indeed, the shear viscosity is one of the coupling coefficients of the effective hydro- 
dynamic description of the theory at lowest momenta and frequency, i.e., , for 

— * 

IR : cj ^ min(T, /i, ■ ■ ■) , |fc| ^ min(T, /x, ■ ■ ■) , (1-14) 

where ■ ■ ■ stand for any microscopic scales of the plasma, other than temperature and 
chemical potential(s) for the conserved charge(s). On the contrary, the microcausality 
of the theory is determined by the propagation of the modes in exactly the opposite 
regime, i.e., , for 

— * 

UV : cj ^ max(T, yU, ■ • •) , |A;| ^ max(T, /i, ■ ■ ■) . (1-15) 

A link between the features of the theory governing its microcausality and its shear 
viscosity is only possible if the same phase of the theory extends over the entire range 



^See [28-30] for further studies. 

''Further work exploring and generahzing this hnk appeared in [32-36] 



of the energy scales — from the infrared to the uhraviolet. In other words, there must 
not be any phase transitions in the plasma. This is precisely what is happening in the 
GB plasma! Since the GB plasma is conformal, and temperature is the only available 
scale in the model, there can not be any phase transition in the theory as a function 
of temperature. The only free parameter of the model is the GB coupling constant 
Xgbi which determines both the shear viscosity ratio and its microcausality properties. 
Hence the link between the two, originally found in [1], is not surprising — rather, in 
a sense, it is an accident. 

Consider a conformal plasma in the presence of chemical potential, defined as a 
holographic dual to appropriately — see the next section for details — generalized GB 
gravity. Assume that this plasma undergoes a second order phase transition below 
some critical temperature Tc oc fi associated with the spontaneous breaking of some 
global U{1) symmetry and the generation of a condensate of some irrelevant operator 

a: 

[=0, T>Tc 

(a) <^ (1.16) 

Clearly, if the model is engineered in such a way that the effective GB coupling of the 
dual gravitational description is 

effective 



A 



GB 



oc a, (1-17) 



it is natural to expect given (1.16) that 



A, 



GB 




uv 

(1.18) 
IR. 



In such a model, microcausality features - governed by the unbroken phase - would 
be completely decoupled from the physics that determines the shear viscosity of the 
symmetry-broken phase of the plasma. Also, as a result, the shear viscosity to entropy 
density ratio in the UV would differ from that in the IR. Thus, although rj/s does not 
flow in any Wilsonian sense (see e.g. [37]), in this construction the decoupling of the 
UV physics from the IR is reflected in the different behavior of rj/s in the two regimes. 
In the next section we present a detailed holographic model of AdS/CFT corre- 
spondence implementing the "decoupling idea" outlined above. We study the thermo- 
dynamics of the model in section 3. The results of the ratio of the shear viscosity to 
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the entropy density and the causahty analysis of the model are discussed in sections 4 
and 5, respectively. We conclude in section 6. 

2 The holographic model 

Following the general idea presented in the introduction, we would like to engineer a 
holographic model of AdS/CFT with a spontaneous symmetry breaking in the IR and 
non-universal shear viscosity in the symmetry-broken phase. 

Our starting point is the holographic model of superfluidity proposed in [38] 
(GHPT) and described by 

/o r \ 3 1 

(2.1) 



r 

'~' super fluid 


r2 /'2T\^ 1 

P T? PM^^ 1 1 Xfiuap p p A \ fSUGRA 
-n. -T^jy-T 1 e -TA^-Tjyo-^p -r ''-scalar ) 

i V d / 4 


with 


rSUGRA _ _1 
^scalar cy 


{d^<f)f + fiin\i^<p{d^e 2A^f cosh^'^ih cosh0) 



(2.2) 

Here, and Q are the modulus and the phase of a complex scalar \1', which is dual to a 
chiral primary operator O with scaling dimension A. By a f/(l) gauge transformation 
we can set ^ = 0. The model (2.1) is a consistent truncation of type JIB supergravity, 
and represents a string holographic realization of the mean-field second-order phase 
transition. 

We now briefly summarize the features and the dynamics of the model (2.1)^. The 
gauge field A^ is dual to a global f/(l) R-symmetry current, and has been normalized 
in such a way that chiral primaries have R-charge \R\ = 2 A/3 [39]. By expanding the 
scalar potential to quadratic order in 0, 

V^(0) = -i - ^^' + • • • ' (2-3) 

we can read off the mass of the scalar rn^L?' = A(A — 4) = —3, and extract the 
dimension of the dual operator. Thus, we can identify \1' with a chiral primary operator 
O of dimension A = 3 and R-charge R = 2. Since the non-normalizable component of 
is set to zero, the dual QFT is a conformal gauge theory. 

Consider this gauge theory at finite temperature T and nonzero chemical potential 
fi. It was found in [38] that for T < T^. ^ 0.0607;U the GHPT plasma undergoes a mean- 
field second-order phase transition associated with the development of the condensate 



^See [38] for further details. 



^Zgf'' ^ C^^alar = "^ [(9,0)^ + A<^^A,A^] + ^ + -^0^ (2.5) 



for O: 

(O) I (2.4) 

[ociT-T,y/\ T<T,. 

On the gravity side, while at high temperatures the background is that of an electrically 

charged AdS black hole, once the temperature drops below Tc the black hole develops 

scalar hair. Since O is charged under the global f/(l) symmetry, the condensation 

breaks this symmetry spontaneously. While the precise value of the critical temperature 

is sensitive to the details of the full scalar Lagrangian C^caUrr"^ (dual to O) in the 

gravitational description, the existence of the transition itself depends only on the set 

of values (-R, A) of the operator in question [40]. Thus, for the purpose of engineering 

the phase transition only — we are going to give up the string theory embedding 

anyway — we simplify 

i[(9,0)^ + 40%A-]+g + A, 

while maintaining (i? = 2, A = 3) for the dual operator O. 

So far, while the simplified Csuperfiuid describes a second-order phase transition, 
the universality theorem of [8] guarantees that the shear viscosity to entropy density 
ratios of the low-temperature (symmetry broken) and the high-temperature (symmetry 
unbroken) phases in this plasma are the same as in (1.1). Furthermore, since in the UV 
the asymptotic geometry described by Csuperfiuid is the same as that of the Reissner- 
Nordstrom black hole in AdS^ [38], the causality properties of the dual plasma must be 
identical to those of A/" = 4 SYM plasma. In particular, we do not expect any violation 
of microcausality^. 

To proceed, we need to introduce higher-derivative gravitational corrections into 
Csuperfiuid in such a way that: 

■ the resulting equations of motion for the background and the fiuctuations are always 
of second order; 

■ these corrections must vanish in the symmetric phase, while being nonzero in the 
symmetry-broken phase; 

■ the phase transition itself should not be destroyed by these corrections. 

A natural modification, obviously satisfying the constraints above, is achieved by 
generalizing the Gauss-Bonnet coupling Xqb of the higher- derivative term in (1.6) as 

^ -^ /3 (^^*)" = /3 02" ^ ^>2, (2.6) 



We explicitly verify this in section 5. 



for some fixed coupling constant /3. Indeed, the Gauss-Bonnet combination leads to 
second-order equations of motion. In the symmetric phase the expectation value of 
Oc = (\1'\E'*)" vanishes, suggesting that the UV properties of the theory must be exactly 
as for (3 = 0; the n >2 condition guarantees that the mass of \E' (and thus the dimension 
of the operator O) will not change as /3 7^ 0. Finally, the sign of (3 will control whether 
the shear viscosity ratio in the symmetry-broken phase is above or below the universal 
result (1.1). 

We can now present our model: 

C = R- —F^,F>^^ + ( — J -t^^^^'^PF^^F.^A, + C,,aiar + Cgb , (2.7) 

where Cscaiar is given by (2.5), and 



^M^pA \ _ 



Cgb = PrL' IR' - AR^^R^" + R^^pxR'^"'^ j . (2.8) 

Thus, while in the UV the scalar field is turned off and one has the simple Einstein- 
Maxwell two- derivative theory, at low energies the scalar field condenses, and controls 
the strength of the higher- derivative GB correction. Note that we have set n = 2 in 
(2.6). Finally, we will consider the dynamics of (2.7) while taking the non-normalizable 
component of to be zero. Thus, C defines holographically a dual conformal gauge 
theory plasma with a global U{1) symmetry. 

To describe the equilibrium state of the plasma - dual to (2.7) - at finite tempera- 
ture and in the presence of a U{1) chemical potential we take the following ansatz for 
the background fields: 

dsl = -cldt^ + cldx"^ + cldr^ , A^ = AS^, 

zovH' zo g . (2-9) 

where {f , g , a , (f)} are functions of the radial coordinate r only. Without loss of gener- 
ality we can choose this radial coordinate such that r = corresponds to the boundary 
while r = 1 is the location of the horizon, i.e., 

lim Ci = , lim C2 = finite = Zo , lim — = 1 , (2-10) 

1 — >-l- r— >-l_ r— >-0+ C2 

for an arbitrary constant Zq. Lastly, we set L = 1. 



(3.1) 



3 The background geometry 

In this section we discuss the thermodynamics of the holographic model (2.7)-(2.9). It is 
straightforward to derive the equations of motion for the background fields {/, g^ a, 0} 
— in the parametrization (2.9), we find two second order equations for {a, 0}, and two 
first order equations for {f,gY^. 

The asymptotic solution near the boundary is given by 

a =q;o + «! r + -p\aQ r^ + 0{r'^) , 

/ =1 + /2 r^ + la\ r' + 0{r') , 
9 =1 - \pI r' + Oir^) . 
Thus, at the boundary the metric reduces to the simple form: 

ds^ = ^{-dt^ + dx^) + ^. (3.2) 

Note that we have set the non-normalizable component of (j) to zero, since we are 
discussing spontaneous symmetry breaking in conformal gauge theories. Altogether 
the UV asymptotics are determined by 4 parameters: {o:o,ai,pi, f2}- Of these, the 
first one, namely a^, is the coefficient of the non-normalizable mode related to the 
U{l)ii chemical potential fi, while the rest are related to the expectation values of 
various operators. The parameter ai, for instance, is the charge density conjugate to 
the chemical potential. 

The asymptotic solution near the horizon y = 1 — r is given by 

a=a1y + 0{y^) , 

^-vUOiy), 

f =0{y) , 

g =g'o + 0{y) , 

where we indicated only the independent parameters. Thus, in the IR altogether we 
have 3 independent parameters: {aJ',Po,5'o}- 



^^These equations are too long to be presented here. They are available from the authors upon 
request. 
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The temperature T and chemical potential fi are 



727r^^ 
/i = 2^0 ClQ- (3-5) 



The thermodynamic potentials are given by 



11 



1 4z^,{aoai - 8/2) 1 4 3 

^^=211 — 3 — ' ^=-Mr°"^ 



(3.6) 



tp ^ ^i.p 

with Q denoting the Gibbs free energy. When translating to gauge theory variables, 

we identify 

1 N^ 

2^ = 8^^' ^^-^^ 

as in the case of A/" = 4 SYM. 

Note that the expression for the entropy density in (3.6) was derived imposing the 
basic thermodynamic relation 

Q = S -sT- fi p. (3.8) 

Alternatively, the entropy density can be computed using Wald's entropy formula [45], 

S = -2n [ d^'-^xV^ -^ e,,ep, , (3.9) 

where E denotes the horizon cross section, h is the induced metric on it and e^,^ is the 
binormal to the horizon cross section. For our geometry the binormal is ttr = C1C3, 
obeying e^^^e^^ = —2, and 

A/* 

2ll — = gf'Pg"^ + 2/30^2 (Rg^Pg'^'^ _ Ag^PR"'' + i?^*^^") . (3.10) 

Putting all the various ingredients together, we find that 



S = -2nAh— e^,ep, = A^ [Atv - 487r/30U 



\2 



(3.11) 

horizon 



where Ah denotes f-^d^ ^x\/—h. Given (2.9), from (3.11) we find that the entropy 
density is 

^\wal,= i3/^4- (3.12) 



^^ These expressions can be obtained following the same procedure as in [43] 
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We mention in passing that a highly nontrivial consistency check on our numerical 
data would be the agreement of the entropy density in (3.6) with the one in (3.12). We 
will return to this point later in this section. 

3.1 Symmetric phase 

In the symmetric phase the field is identically zero, which tells us that the parameters 
Pi and Pq vanish. Thus, once {T, fi} are fixed — i.e. given {zq, ag} — "we are left with 
4 integration constants {ai, f2,ai,gQ} — precisely the correct number necessary to 
uniquely solve a coupled system of 1 second-order differential equation (for a) and 
2 first-order differential equations (for {f,g}). Actually, in this case the background 
equations of motion can be solved analytically^^. We find: 

a = ao{l-r), / = l _ r^ + ^(r^ - r^) , g = l. (3.13) 

In this case the thermodynamics is that of the A/" = 4 SYM plasma with the same 
chemical potentials for all the U{1)^ C SU{4) R-symmetry global charges [46]. The 
expression for the temperature reduces to 

T.5(:-f), ,3.14) 

in terms of which the entropy density becomes 

' '-'^' . (3.15) 



211 



3.2 Broken phase 

In the broken phase the field is no longer zero, and the parameters Pi,Pq are now 
turned on. Thus, for a fixed {T,fi} — given {zo,ao} — we are left with 6 integra- 
tion constants {q:i,Pi, f2,o.i,Po^go}- This is precisely the correct number necessary to 
uniquely solve a coupled system of 2 second-order differential equations (for {a, 0}) 
and 2 first-order differential equations (for {f,g}). We use numerical techniques de- 
veloped in [44] to study the thermodynamics of the low-temperature symmetry-broken 
phase of (2.7) for different values of the coupling constant f3. 
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Our numerical results are in excellent agreement with the exact analytical result. 



12 



8 Q 



^blue/^ 



symmetric 




Figure 1: (Colour online) Left plot: the Gibbs free energy densities of the symmetric 
phase (purple curve) and the symmetry-broken phases at /3 = (black curve) and 
/3 = — 1 (blue curve) as a function of ^^. Right plot: the ratio of the Gibbs free 
energies in the symmetric and the broken phases at /3 = — 1. 



We find that the mean field second-order phase transition at /3 = persists for 
/3 7^ 0. The positive values of (3 tend to increase the ratio - (for fixed non-normalizable 
modes {2:0, ao}), while the negative values of (3 tend to decrease it. We have not per- 
formed an exhaustive analysis of the phase diagram of the system, but rather identified 
interesting values of the coupling with regards to the ratio of shear viscosity to entropy 
density (see section 4). 

A representative case of this analysis is the comparison between the Gibbs free 
energy densities of the broken and unbroken phases as a function of ^^ at /3 = {0,-1}, 
which is shown in Fig. 1. On the left, the "thin purple" curve represents the Gibbs 
free energy density of the symmetric phase 

1024 3x2 + 4 - xv/9x^Tl6 



8 Vt 



7]-2jY2 y4 — J purple 



fp 



27rT 



X 



(3.16) 



/u y 27 x^{V9x^ + 16 - 3x)^ 

On the other hand the "thick black" and the "thick blue" curves represent the Gibbs 
free energy densities of the broken phase at, respectively, (3 = and (3 = —1. As the 
temperature increases, the condensate of the dimension-3 operator (O3) dual to the 
holographic scalar in (2.5) 

{0^)^Vi (3.17) 

decreases, ultimately vanishing at some critical temperature Tc = Tc{(3). We find 



27rT, 



0.65396(3) 



13=0 






0.63040(9) 



(3.18) 
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Figure 2: (Colour online) Ratios of the Gibbs free energy densities in the broken and 
unbroken phases as a function of ^^, for select values of the coupling. From top to 
bottom: (3 = —1 (blue), /3 = — 2 (red), /3 = — 5 (orange) and /3 = —10 (green). 



The right plot in Fig. 1 represents the ratio of the free energies in the broken and un- 
broken phases at /3 = — 1. Clearly, the broken phases are thermodynamically favorable 
at low temperatures. 

To get to the interesting regime in the shear viscosity ratios - (in the broken 
phases) we need to get to temperatures several times smaller than the appropriate 
critical temperature. Fig. 2 shows ratios of the Gibbs free energy densities in the 
broken and unbroken phases for a select set of couplings, /3 = { — 1, —2, —5, —10}. The 
broken phases, while being closer and closer to the unbroken phase as /3 decreases, are 
thermodynamically preferable for each given temperature. 

An important consistency check on our numerical analysis is the comparison^^ of 
the entropy density derived from the basic thermodynamic relation (3.8) — see (3.6) 
— and the one obtained directly from the Wald's entropy — see (3.12). In all instances 
we find 

(3.19) 



S\Wald 



1 



< 10 



-6 



13 



Of course, this is a nontrivial check in the symmetry broken phases only. 
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4 Shear Viscosity 

In the hydrodynamic approximation, linear response theory imphes that the retarded 
Green's function of the stress energy tensor of the conformal fluid in the tensor channel 
is given by [47] 

G™ (tt) = ^,^ = O) = -tjdtdxe'^'eit) {%yix),T,yiO)]) 

= p(l- 2.tt,. ^- Z±lpJ^ + o(^^ ' ^"''^ 

\ s P 

Techniques for computing this correlation function in a dual gravitational model are 
well developed [12], and we will not review them here. In particular, the analysis which 
we perform in this note is equivalent to that in [16]. 

On the gravity side, computation of the Green's function (4.1) entails adding a 
metric perturbation of the form 

yxy ' yxy ~r it-xy ; \^'^) 

and finding the effective action for the fiuctuation h^y. Thus, we take the metric to be 
dsl = -cldf + cl{dx^ + dy'^ + dz'^ + 2e^ dxdy) + cldr"^ , (4.3) 

and expand 

*(t,r,z) = y"^e-*+'^>,(r). (4.4) 

Since we are ultimately interested in the correlator at vanishing spatial momentum, 
we can set A; = at this stage and consider perturbations which depend on (r, t) only. 
Expanding the action (2.7) to second order in the perturbations, we can easily see that 
the effective action for the fluctuation is of the form originally found in [16], 



c(2) = 1 ! ^dr 



^ 21% J (27r) 



A{rWLV-k + B{r)^'^ip'_^ + C{r)^'kip_k + D{r)ipkip-k + 



+E{r)^lip"_k + F(r)v9^V-fcJ + ^GB + lC,ounter ■ (4.5) 

Here ICgb denotes the generalized Gibbons-Hawking boundary term, needed to en- 
sure a well-deflned variational principle, and ICcounter is a local boundary counterterm, 
necessary to remove UV divergences in the stress energy tensor correlation functions^^. 



^^For the construction of the generahzed Gibbons-Hawking boundary terms and counterterms for 
Einstein-Maxweh theory in the presence of generic R^ corrections, see [51]. 
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For our model (2.7) we find that E{r) = 0, which is in agreement with the expec- 
tation that the equation of motion corresponding to Gauss-Bonnet gravity should not 
contain more than two derivatives. When E = the generalized Gibbons-Hawking 
term takes the simple form 



/C 



GB 






-^vW-k - -^vW-k 



(4.6) 



Furthermore, the local boundary counterm is precisely as in the case of pure AdS^ [48] : 



/C, 



counter 



dM 



d^k 

(2^ 



B 



(4.7) 



It turns out to be particularly convenient to rewrite the action in the form 



q(2) _ ^ 



2^p JM 



dr 



B-A-^)vW- 



9/3 , 

^'p JdM 



d^k 



n C" - A" 

^ z ]<Pk'^-k 



+ 



(4.^ 



from which one can easily read off the radial canonical momentum for the scalar ip, 



UJr] 



6S, 



(2) 



1 / p'^ 

^(B-^-yK. 



Introducing an "effective mass" term for the scalar fluctuation 

C'-A"^ 



M 



e//l 



k'^ 



the scalar equation of motion can be written in the simple form 

drilk = Meffipk ■ 



(4.9) 



(4.10) 



(4.11) 



By making use of the background equations of motion it is straightforward to verify 
that Meff = (9(ti3^), which in turn means that the radial flow of 11 is trivial in the 
tt) — )■ limit, and the mass term M^ff does not contribute to (4.1) to order (9(lt)). 

Finally, we note that evaluating the on-shell action to order 0(tt>) turns out to be 
equivalent to evaluating the following boundary term^^ 

d^k 



5". 



on-shcU 



2/3, 



d^k 



(2vr)' 



B 



J^k 



A - ^F'] ^i^_k + lie -A' + B)vkV-k 



r=l 



r=0 



,(4.12) 



^^Explicit expressions for {A, B, C, F, B} arc too long to be presented here. They are available from 
the authors upon request. 
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with the flux J-^ directly related to the retarded Green's function: 

Gf,,,, = - hm .Y' ( ^ ■ ^4-13) 

Much like the background in the broken phase (see section 3), the equation of 
motion for the scalar ipk (4.11) has to be solved numerically. The fluctuation ipk must 
satisfy an incoming wave boundary condition at the horizon [49] 

ipk = {l- r)-^"/2 Mr) , lim^fc = 1 , (4.14) 

r— ^1 

where ipk{r) is regular near the horizon, r — ?> 1_. Note that we used a conventional 
normalization for ipk. To compute the correlator (4.13) to order 0{w), we need to solve 
(4.11) to order O{vo) as well. We represent 

Mr) = i^lir) + itt) Mr) + C(m') . (4.15) 

Demanding regularity at the horizon we find 

^fc = 1 (4.16) 

identically. The second order linear inhomogeneous equation for ipl has the following 
asymptotic solution 

^l =^0 - ^^ + ^2 r^ - ^r^ - Q + 1 /2 + ^ h^2) r' + 0{r') , 

^l =0{y) . 

close to the boundary r — )■ and the horizon y — )> correspondingly. It is uniquely 
specified by two parameters {iIjq,iIj2}- Comparing the holographic expression for the 
Green's function (4.13) with that of the hydrodynamics (4.1), we arrive at a fairly 
simple expression for the ratio of shear viscosity to entropy density: 

T] 3 1 + 4^-2 



s Stt aotti — 3/2 
4.1 Shear viscosity of the symmetric phase 



(4.18) 



In the symmetric phase the background is known analytically, and is given by (3.13). 
The equation of motion for ipl takes the following form 



4a2(r3_r2)-9(l + r2) 2alr^ + 9 

{r^ - r){2ay - 9{r + 1)) ^^ ~ 2{r^ - r){2ay -9{1 + r)) 
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Figure 3: (Colour online) Ratios of shear viscosity to entropy density in the broken 
phase for select values of the coupling: /3 = — 1 (blue), /3 = — 2 (red), /3 = — 5 (orange) 
and /3 = — 10 (green), as a function of ^^. The dashed black line indicates the Gauss- 
Bonnet viscosity bound: rj/s > 16/25 [1]. 



which can be solved analytically: 



V', 



dx 



lalx + 9 



2(9(1 + x) -lalx'^) 



From (4.20) we can extract 



V'2 



1 

4 



an 



(4.20) 



(4.21) 



Finally, using the explicit solution (3.13), relation (4.18) reproduces the shear viscosity 
of the A/" = 4 SYM plasma in the presence of t/(l)/j chemical potential [50] : 



s 



1 

47r 



(4.22) 



4.2 Shear viscosity of the broken phase 



In the broken phase the equation for motion for tl}\ (4.15) and the shear viscosity to 
entropy density ratio (4.18) must be computed numerically. First, for /3 = we find 



4vr!^-l 



< 10 



-7 



(4.23) 



/3=0 



and therefore recover the universal t]/ s = 1/4tt result [8] expected for a two-derivative 
theory. 

Fig. 3 represents the results of the numerical analysis of the shear viscosity in the 
symmetry-broken phase of the holographic model (2.7), for select values of the coupling. 
We show /3 = — 1 (blue), /3 = — 2 (red), /3 = — 5 (orange) and /3 = — 10 (green), as a 
function of ^^. Notice that for /3 = — 1 the shear viscosity remains above the causality 
bound for Gauss-Bonnet gravity found in [1] 

Tl 16 
47r-> — = 0.64, (4.24) 

s 25 

while for the other values of /3 we consider, it dips below this bound. As we observed 
in the thermodynamic analysis of the broken phase, see Fig. 2, decreasing /3 makes 
the broken phase (while still thermodynamically preferable) closer and closer to the 
unbroken phase. Correspondingly, for smaller values of /3 the shear viscosity is closer 
to the universal result down to lower and lower temperatures. However, for sufficiently 
low temperatures it drops even steeper. It is technically challenging to perform our 
numerical analysis reliably at temperatures lower than those reported; nonetheless, 
the data obtained suggests that the holographic plasma (2.7) does not have any lower 
bound on the ratio of the shear viscosity to the entropy ratio as one varies f3. 

Finally, we note that for positive values of (3 the broken phase of the holographic 
plasma (2.7) has a shear viscosity ratio exceeding the universal result; we have not 
studied this parameter regime in detail. 

5 Causality of holographic superfluid plasma 

Consistency of a holographic plasma as a relativistic quantum field theory requires that 
it does not propagate modes faster than the speed of light. The dispersion relation 
of the linearized fiuctuations in the plasma is identified with the dispersion relation of 
the quasi-normal modes of a black hole in the dual gravitational background. There 
are three types of quasi-normal modes in gravitational geometries with translationally 
invariant horizons [52]: 

■ a scalar channel (helicity-two graviton polarizations); 

■ a shear channel (helicity-one graviton polarizations); 

■ a sound channel (helicity-zero graviton polarizations). 

In the case of the GB plasma, the lower bound on the shear viscosity (the upper bound 
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on Xgb (1-11)) comes from the scalar channel quasi- normal modes [1]. On the other 
hand the upper bound on the shear viscosity (the lower bound on Xgb (1-12)) comes 
from the sound channel quasi- normal modes [27]. In our case, the study of the quasi- 
normal modes in the sound channel is the most difficult — it requires understanding 
holographic viscous hydrodynamics in the presence of Goldstone modes associated with 
the spontaneous breaking of a global f/(l) symmetries. To our knowledge such theory 
has not been developed yet^^. Instead, as in [1], we limit our discussion to the scalar 
channel quasi-normal modes. We expect that analysis of the other channels will not 
change our conclusions with regards to causality. 

Our discussion here follows closely [1]. Due to their complexity, we omit most 
technical details^^. The quasi- normal equation for the scalar channel fluctuations for 
the holographic plasma dual to (2.7) takes the form 

^[sca«ar](^) + ^scalar ^[sca«ar](^) + ^scalar Z[scalar]{r) = . (5.1) 

Following [1], it is possible to introduce a new radial coordinate y = y{r), with y — )■ — oo 
corresponding to the horizon and y — )> 0_ corresponding to the boundary, and to rescale 
the radial profile as 

■^[scalar] ^ "^ "^[scalar] ; V"-*'^/ 

such that (5.1) can be brought into the form of an effective Schrodinger equation: 

y ^[scalar] ~r tV [scalar] "V^ [scalar] C^ "V^ [scalar] ; TT ' ^^ "iT 

where f/[scalar] = f^fscalar] + ^ ^^[scalar] > ^c^ f^[scalar] = 5 "^c, f/[scalar] 7^ • 

(5.3) 

Notice that in the limit fc — ?> oo (or /i — )■ 0), everywhere except in the tiny region 
y > — ^ the dominant contribution to the effective potential U scalar comes from U^^aiar- 
Thus, in this limit it is a good enough approximation to take 



^' f/[scalar] = <^ ' (5-4) 

l+oo y>0. 

As explained in [1,23], the bound states of the resulting 1-dimensional quantum 
mechanical problem (5.3) with "energy" Cg > 1 point to the presence of quasi-normal 



^^For a step in this direction see [53]. 

""^^Oniitted expressions are available from the authors upon request. 
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Ui 



[scalar] 




Figure 4: (Colour online) Effective scalar potential in the symmetric phase, see (5.5). 
The blue, green and red curves correspond to, respectively, ao = {1, 3, 5}. The regime 
ao > 3 is unphysical, since it corresponds to negative temperatures. 

modes in the plasma, propagating faster than the speed of light. On the other hand, 
bound states with energy c^ < 1 indicate the presence of instabilities (tachyonic modes 
in the plasma in the limit A; — )■ oo). 

5.1 Causality of the symmetric phase 

The effective potential U[scaiar] defined in (5.3) and (5.4) can be computed analytically 
in the symmetric phase, and is given by: 



Ui 



[scalar] 




(1 + lal) + |«2 
r < . 



<r < 1 



(5.5) 



Fig. 4 shows the potential U[scaiar] above for select values of ao- The blue, green and 
red curves correspond, respectively, to ao = {1,3,5}. In the ao = 1,3 plots, U[scaiar] 
decreases monotonically between the boundary and the horizon, and never develops 
a maximum in the intermediate region. In all cases there are no bound states with 
energy c^ > 1 - there are no superluminal quasi-normal modes and, as expected, the 
theory is causal. Notice that for ao > 3 the potential (5.5) develops a negative energy 
minimum, which implies the existence of negative energy bound states, and as a result 
tachyonic (unstable) quasi- normal modes in the fc — ;■ oo limit. However, this does not 
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Figure 5: (Colour online) Effective potential U[scaiar] in the symmetry broken phase of 
the holographic plasma (2.7) at /3 = —5 and ^^ = 0.01(9), corresponding to An- = 
0.4(9) . The dashed red lines correspond to c^ = {0, 1}. 

cause any problems: for a^ > 3 the temperature of our plasma in the symmetric phase 



T 






1 - 



Or 



(5.6) 



becomes negative, and the tachyons are therefore not physical. 

5.2 Causality of the broken phase 

In the symmetry-broken phase the effective scalar potential U[scaiar] (5.3), (5.4) can only 
be computed numerically. A representative example of such computation is shown in 
Fig. 5. Here, /3 = -5 and ^^ = 0.01(9), corresponding to Air^ = 0.4(9) (the low 
temperature endpoint of the orange curve in Fig. 3). Notice that this potential does 
not support bound states with energy cj. > 1; neither does it support states with c^ < 0. 
We conclude that, at least in the scalar channel, the gauge theory plasma holo- 
graphically dual to the gravitational model (2.7) does not violate causality. We also 
find that it does not contain any tachyonic modes in the fc — ^ oo limit - the theory 
appears to be perfectly well-behaved over the entire range of parameters. Unlike the 
case of the GB plasma (1.6), here the self-consistency of the CFT doesn't place any 
constraints (whether from below or above) on the size of the higher-derivative cou- 
pling. While it is possible that the shear and sound channels might lead to additional 
instabilities, previous studies [27, 34] suggest that this should not be the case. 
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6 Conclusion 

In this paper we have argued that microscopic constraints (causahty, positivity of 

energy, etc.), while important for the general consistency of a plasma as a relativistic 

quantum field theory, are not necessarily responsible for setting the lower bound on 

the ratio of shear viscosity to entropy density in the plasma. The basic reason is that 

the hydrodynamic transport of the system is determined by its infrared properties, 

which do not necessarily enter into the microcausality analysis of the theory. To this 

end, we generalized the holographic model of "GB plasma" introduced in [1] in such 

a way that the Gauss-Bonnet coupling of the former is replaced with an (irrelevant) 

operator. Our holographic model, see (2.7), undergoes a second order phase transition 

at low temperatures, where this operator develops a vacuum expectation value. As 

a result, the effective Gauss-Bonnet coupling in our model is nonzero in the broken 

phase (which is necessary to generate the non-universal ratio of shear viscosity to 

entropy density), but being identified with an irrelevant operator it does not effect the 

ultraviolet properties of the model — the dynamics at high energies is equivalent to 

that of holographic superconductors [38]. 

We identified parameters in our model where the shear viscosity drops below the 

causality bound^^ 

r? 1 16 , , 

7^4^ 25- <«^1' 

determined in [1]. It would certainly be interesting to identify the lowest bound in our 

model — however, this is not the main focus of this paper. It is clear that, whatever 

the lowest bound (assuming it exists) on the shear viscosity ratio in holographic plasma 

(2.7), it does not affect its causal properties. To complete the analysis one would need 

to study causality in the vector and the sound channels of the plasma quasi-normal 

spectrum-*^^ [27]. As we already stated, we do not believe that such analysis would 

modify the physical picture presented here. 

To summarize, the question of the bound on the ratio - suggested by a quasi-particle 

picture of the fluid, its very existence, and the physics that determines it remains open. 



"'^^The current lowest bound on the ratio of the shear viscosity to the entropy density in 4-diniensional 
plasma was reported in [54] . It is not clear though whether the model discussed there is a consistent 
relativistic QFT. 

^^The analysis of the sound quasinormal models are most challenging and will be reported elsewhere. 
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